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INTERACTIONS BETWEEN EHRENFEST’S URNS ARISING FROM 

GROUP ACTIONS 

HIROSHI MIZUKAWA 


Abstract. Ehrenfest’s diffusion model is a well-known classical physical model consist¬ 
ing of two urns and n balls. A group theoretical interpretation of the model by using the 
Gelfand pair (Z/2Z;S'„, Sn) is provided by Diaconis-Shahshahani [5]. This interpretation 
remains valid for an r-urns generalization, in which case, the corresponding Gelfand pair 
is {Sr I Sn,Sr-i I Sn). In these models, there are no restrictions for ball movements, 
i.e., each balls can freely move to any urn. This paper introduces interactions between 
urns arising from actions of finite groups. It gives a framework of analysis of urn models 
whose interaction between urns is given by group actions. The degree of freedom of 
ball movements is restricted by finite groups actions. Furthermore, for some cases, the 
existence of the cut-off phenomenons is shown. 


1. Introduction 

Ehrenfest’s diffusion model [3j is a stochastic process consisting of n-balls {6i,..., bn} 
and two urns {Uq, Ui}. The model is constructed in the following way: Let 'L/2'L = {0,1} 
be a cyclic group of degree 2. We set U{2,n) = {(xi,... ,x„) | Xi G Z/2Z}. We identify 
an element (xi,..., x„) of U (2, n) with a state indicating whether the ball 6* is in the urn 
Uxi {xi = 0 or 1). We dehne a stochastic matrix on U{2,n) by 

{ l/(?7, + l) x = y, 
l/{n + l) d{x,y) = l, 

0 otherwise. 

Here, d{x,y) = \{i \ Xi ^ yi, 1 < i < n}\ is the Hamming distance on U{2,n). The second 
condition above means that each ball moves to another urn under a probability 1 in n + 1 
at each step. We can consider a generalization of the above-mentioned stochastic process 
above by increasing the number of urns. Let U{r,n) = {(xi,...,x„) | Xj G 
Here, Sr is the symmetric group on {0,1, • • • ,r — l} and Sr-i = {a G S'r | a{r — l) = r —1} 
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is a subgroup of Sr- Then, we define a stochastic matrix on U{r,n) by 

f iiTT ^ = 

= { (r-l)(n+l) = 

[O otherwise. 

Note that setting r = 2 gives the original 2-urn case is recovered. In this generalization, 
note that each ball move to any urn in each step. 

Diaconis and Shahshahani [2] analyzed the asymptotic behavior of the iV-step prob¬ 
ability of the original Ehrenfest’s diffusion model. In particular, they showed that the 
cut-off phenomenon occurs in this case. For general r-urns cases, Hora mm gave precise 
analysis of the asymptotic behavior of the iV-step probability. He also showed that the 
cut-off phenomenon occur in these cases. 

In these cases, each model is realized by certain finite homogenous space arising from 
a hnite Gelfand pair {Sr I Sn, Sr-i I Sn) {n = 2, 3,4,...). We remark that the theory of 
hnite Gelfand pair works well in these cases. One can find detailed descriptions of the 
theory of Gelfand pairs in Macdonald’s book [9] and its relation to probability theory in 
[1]. However, in these cases, balls can freely move to any urns at each step, or in other 
words, there are no interactions between urns. Of course, many other interactions could 
be considered. In this paper, we give some interactions between urns through actions of 
finite groups actions. It gives a framework of analysis of urn models whose interaction 
between urns is given by group actions. 

The main tool considered in this study is a multivariate hypergeometric type orthogonal 
polynomial called a multivariate Krawtchouk polynomial [U [H]. An important feature 
is that the Krawtchouk polynomials can express the zonal spherical functions of a finite 
Gelfand pair {K lSn,LlSn) [ID] . Multivariate Krawtchouk polynomials are known to have 
an application in the stochastic model of the poker dice game [5] . In this paper, we present 
another stochastic theoretical application of the multivariate Krawtchouk polynomials, 
that is, to the Ehrenfest’s diffusion models. 

The remainder of this paper is organized as follows. In section 2, we construct interac¬ 
tion urn models by using a finite Gelfand pair {K, L) and prepare some tools to analyze 
our models. In section 3, we compute upper and lower bounds of the total variation 
distance between an Wstep probability and the uniform distribution on our models. As 
a conclusion of our discussion, we show that our models have the cut-off phenomenon. 

2. Preliminaries 

Let iL be a finite group and L be its subgroup. Set = |^ ^ which is an element 

of the group algebra CK. Then the following conditions are eqivalent: (1) The Hecke 
algebra 'H{K,L) = e^CKei C CK is commutative. (2) The permutation representation 
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CKei is multiplicity-free as ii'-module. These condition are satisfied, the pair {K,L) 
is called a Gelfand pair. Here we remark that 'H{K, L) can be identihed with the ring 
of bi-L invariant functions, i.e. functions on double coset L\K/L, on K. Assume from 
now that {K, L) is a Gelfand pair and multiplicity-free decomposition of the permutation 
representation is given by Ce^ = 0*rQ 14. It is fact that s = dim?^(A', L). Let Xi be 
the character of I 4 . We dehne function oji by uji{g) = The functions 

cuj (0 < z < s — 1) are called zonal spherical functions of Gelfand pair {K, L). The zonal 
spherical functions is a basis of 'H(iL, L) satisfying orthogonality relations 

^ = ^ij dimld-b 

I I geK 

Therefore we can expand a bi-L invariant function on K by using the zonal spherical 
functions. One can hnd more detailed descriptions of the theory of Gelfand pairs in 
Macdonald’s book [9] 

Let K be a hnite group. The wreath product KlSn = {(s'!, • • •, 5'^; cr) \ Pi & K, cr E S'„} 
is the semidirect product of iL" with Sn whose product is dehned by 

{gi,.. .,gn;a){g [,... ,5'nW') = • • • > 

Let a positive integer t be the number of conjugacy classes of K. The irreducible repre¬ 
sentations oi K I Sn is determined by the t-tuples of partitions (A°,..., such that 

|A°| -I--h |A*“^| = n (see [8]). 

The next proposition is the fundamental result of BDl. 

Proposition 2.1 ([iQ]). If{K,L) is a Gelfand pair, then {KlSn,LlSn) is also a Gelfand 
pair. The double coset and the zonal spherical functions of {KlSn, LlSn) are parameterized 
by X{s,n) = {(/cq, • • •, kg-i) \ ko + ■ ■ ■ + kg-i = n}. The values of the zonal spherical 
functions Qk at x = (xi,..., Xn] a) E K I Sn are given by the coefficients of t^ in 

\j=0 / kGX{s,n) 

(1) = for k = {ko,..., kr-i) E X{s,n). 

The functions Qk in Proposition 12.II are called the multivariate Krawtchouk polynomials 

(mum)- 

3 . Formulation of interaction urn models arising from group agtion 

Throughout this paper, a pair of hnite groups {K, L) is a Gelfand pair. For the ir¬ 
reducible decomposition of permutation representation CKcl Vi, let Vo be the 

trivial representation of K. Gonsider di = diniFj (0 < i < s — 1). Let cu* be a zonal 
spherical function of {K, L) afforded by V^. 


^{K,L-,n){x) - 
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Fix an element Xq E K — L. Set Dq = LxqL. Put 


r = 


m = 


\K\ 

\L\ > 

_ tPoi 
FI ■ 


For X G IF"', we denote by x a left coset xL^ G JF"'/L"'. 

Definition 3.1. Let 0 < mp < 1. 

(1) We define a function on {K/L)^ x [K/LY by 

1 — mp x~^y G L” 
P{x, y) = { p/n 


0 


X ^y E X Dq X ^ ^ 

otherwise. 


(2) We define a matrix P(K,L-,n) by 

P{K,L-,n) = {p{x,y))x,ye{K/L)"x{K/L)^- 

In the above dehnition, we remark two points. First the conditions x~^y G L” and 
x~^y E X Dq X do not depend on the choice of a complete set of representatives 

of {K/LY. Second the random walk depends on the choice of xq- Here, if we identify 
{K/LY = • • •, Xn) I Xi E K/L} with the states U{r,n) = {(xi,..., Xn) | 0 < Xj < 

r — 1} as a set, we can interpret Dehnition 13.11 as a certain stochastic process of n-balls 
and r-urns. The interaction between urns is determined by a Gelfand pair {K,L). To 
understand the meaning of the interpretation, we provide some examples. 

Example 3.2. In the below examples, we define d(x,y) = |{i | x* Y yi}\ x = 
{xI,...,xL) andy={yI,...,yL). 

(1) Let K = Sr and L = Sr-i- Here, Sr-i is a subgroup of Sr generated by transposi¬ 
tions {{i, i + 1) \ 1 < i < r — 2}. A complete system of representatives for the left 
cosets and that for the double cosets of L in K are given by 

/F/L = {ei^, (l,r),(2,r),...,(r-l,r)} and L\JF/L = {e^^, (1, r)}. 

Let Xq = (1,?"). Then, we have |Do| = (r — l)(r — 1)! and m = {r — 1). We 
Dq i Y j 

Now, we identify a left coset {i,r)L with an urn 

L i = j. 

Ui {0 < i < r — 1). If we set p = , then P{K,L;n) is the same matrix of an 

r-urns generalization of Ehrenfest’s diffusion model introduced in Section\^ 

(2) Let JF = (a I a'' = e) be a cyclic group and its subgroup L = {e}. Setting xq = a, 
we have m = |LaL|/|L| = 1. Then, 


have {i,r){j,r) E 


a E Dq j — i = 1 (mod r). 
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Set X = X = (a*\ ... ,a*") and y = y = Then, the {x,y)-elements of 

P{K,L-,n) are given by 

{ 1 — p X = y 

p/n d{x,y) = 1 and jm — = 0 or 1 (mod r) 

0 otherwise. 

This case can be interpreted as an n-balls and r-urns model as follows: At each 

step, a ball is randomly chosen with the probability 1/n and the chosen ball, which 
is in urn Ui, either moves only to its left neighbor urn Ui^i + 1 G Z/rZ) with 
probability p or stays in the same urn with probability 1 — p. 

(3) We consider a dihedral group K = Dr = {a, h \ a'^ = b^ = {obY = 1) and its sub¬ 
group L = (b). A complete representatives of L\K/L is given by 
Set xq = a. We have m = \LaL\/\L\ = 2. Set x = (a*S...,a*") and y = 

,..., aP). Since G LaL ^ j = i ± 1 (mod r), the {x,y)-elements of 

P{K,L-,n) are given by 

{ 1 — 2p X = y 

p/n d{x,y) = 1 and = 0 or ±1 (mod r) 

0 otherwise. 

This case can be interpreted as an n-balls and r-urns model as follows: At each 

step, a ball is randomly chosen with the probability 1/n and the chosen ball, which 
is in urn Ui, moves only to either of its both neighbor urns Ui-i and Ui+i {i, i ± 
1 G Z/rZ) with probability p or stays in the same urn with probability 1 — 2p. 
(i,i ± 1 G ’L/r'L). 

(4) In general, {K x K, AK) is a Gelfand pair for any finite group K. Here AK = 
{{9^9) I 9 £ K} is a diagonal subgroup of K x K. This fact supports that our 
model includes numerous examples. 

From these examples, we see that m gives the number of directions of ball movements. 
We dehne an action of iF ^ S'„ on (iF/L)"' by 

i.9li • • • 1 9n • (^1) • • • ) ^n) (.9Wcr{l) y • • • t 9rHij{n )) •> 

where xTi = XiL (xj G K). From this dehnition, we see that K I Sn acts on (K/L)^ 
transitively. 


Proposition 3.3. P{K,L-,n) is a K I Sn-invariant stochastic matrix. 
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Proof. It is clear that p{gx, gy) = p(x, |/) > 0 from Definition 13.11 Let ok be the identity 


element of K and let e = (ci^,..., ok) G {K/LY. We compnte 



y&{K/L)' 


y£(K/L}^ 




y£(KIL}^ 


n 


p{e,e) + E E p{e,y) 


y&pKp ^xDo/Lx{eKp * 
(1 — mp) + n{mp/n) = 1. 


□ 


Thus, with the natural identihcation between (K/L)'^ and KiSn/LlSn, we can analyze 
our stochastic model by using the theory of a Gelfand pair {K I Sn, L I Sn). 


4. iV-STEP PROBABILITY 


We attempt to obtain the iV-step probability, i.e., we compute P^ L-n)- recall a 
brief review of the relation between the iV-step probability and the theory of hnite Gelfand 
pairs in our case. We £x e = (Pk, ■ ■ ■ ^Pk) G {K/LY as an intial state. Let p(x,y) be an 
(x,^)-entry of P(K,L-,n)- Set u{x) = p{e,x) and i>{g) = -^f^ig) {g e KlSn), where i = \L\. 
Then, from Proposition 13.31 h is a bi-L I Sn invariant function on K I Sn- Therefore, we 
can expand h as a linear combination of D^s. Let h = Ylikex(sn) \kis where 

d, ^ 4° ■ ■ ■ dt-pi) is the dimension of an irreducible representation that affords to ilk 
([10]). Let k'Nig) be the (e, ^)-entry of P^k L-n)- compute 




XI,---,xm-i&(K/LY 



L I SnY*^{g) (A^th. convolution power). 


Then the idempotence of the zonal spherical functions [9] gives 



k£X{s,n) 


We call the /(A;)’s the Fourier coefficients. 
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Theorem 4.1. For k = {ko,..., kg-i) € X{r, n), we have 


's-l 


ki 


f(k) = (1 - mp) +mpl^ j • 


. i=0 


Proof. We must compute the value of Qk evaluated at x G x Do x ^ Namely, we 
compute the coefficient of in ^(K,L-,n){x) = ^ j • Through 

direct computation, we obtain 


S—1 


h 


^k{x) = —Ui{xo). 


i=0 

Then, the Fourier coefficient is computed as follows: 
/(^) = ^(9)^ki9) 

geKiSn 

n—1 


1 .. . . Ip 

i'^nl n 


i^nl 


(1 — mp)(Fn\ + 


E E 


VLk{g)n\ 


j=0 g£L^ xDoxL”-J-^ 


s—1 


ki 


= (1 — mp) + mp —cuhxo). 

^ n 


1=0 


Here, we use m^ = |Zlo| = |TxoT| at the last equality. 

5. Upper and Lower Bounds Evaluation of 


□ 


5.1. Upper bound evaluation. Let tt be a uniform distribution on (K/L)"'. We at¬ 
tempt to evaluate the total variation distance \ \uj\f — ttHtv for some cases constructed by 
Dehnition 13.11 

Definition 5.1. Set M = max{a;j(xo) | 1 < f < p — !}■ 

For the remainder of this paper, we assume the following three conditions: 


Assumption 5.2. 

(2) M < 1. 

(3) 0 < mp <1/2. 


(1) a;i(xo) G M (1 < i < r — 1). 


Example 5.3. Example\3fE- (1) and Example\3fM-(3) satisfy Assvmption \5. 2 - f 1) and (2). 


Theorem 5.4. Set N = 


2mp{i-M) (fogri(p — 1) I- c) . Under Assumption \ 5. 2i we have 


||pv-vr||Tv< ^ (exp (exp (-c)) - 1 ). 

Here H^at — vtHtv is the total variation distance between un and the uniform distribution 

71 . 
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Proof. We set • • • k'l’Lf for k G Xi(s, n) and d = {do,..., dg-i). We use the 

following facts for our computation below: 

(a) the upper bound lemma [H Corollary 4.9.2]. 

(b) di d2 dg—i = T — 1. 

(c) exp {—x) > I — X, (x < 1). 

Let 


'ip{ko) = max 


's-l 


ki 


1 + mpi^ —{uJi{xo) - 1 ) 


. i=l 


I /ci H- h /Cs-1 = n-ko 


for 0 < /cq < u — 1. Since mp <1/2 and —2 < M — 1 < 0, we have 


■j/(/co) = 1 + mp ^ — 1). 


n 


Let Wi(s, n) = X{s,n) — {(n, 0,..., 0)}. We compute 


\Wn — TtIItV < T '^2 ^ 

(a) 4 


kGXi (s,n) 


kGXi (s,n) 


n 

rP 


S — 1 


(1 — mp) + mp I '^2 —^i{xo) 


.1=0 


2N 


1 + mp 


- 1 ) 


2N 


< 


n—1 

iS: 

ko=0 


kGX(s,n—ko) 


n — k( 


P{ko) 


2N 


n—1 


{b) 4 


ko=0 


E L 


1 ^ 

“ 4 ^ {n — ko)\ 

ko=0 ^ 


(r — 1)"' ^° ( 1 — mp -^(1 — M) 


n — ko 


n 


2N 


n — kn 


(r — 1)"' ^°exp ( —2Nmp -(1 — M) 


n 


1 

< - 

(c) 4 


^ fco ^ 1 /^ '\ 

^^ 2 ^ ~ 1)^° 6xp 2A^mp—(1 — M)J 


Here, we set 


N = 


n 


2mp{l — M) 


(logn(r — 1) + c). 
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This leads to 

T ^ ^0 1/^ 

- -(r — 1)^° exp (—2iVmp—(1 — M)) 

4 knl n 

ko=l ^ 

^ ko 

- 4 ~ 1)*^° exp (-fco(logn(r - 1) + c)) 

fco = l 

_ 1 exp (- 0 )"^° 

“ 4 ^ k\ 

ko=l 

= ^ (exp (exp(-c)) - 1). 


□ 


5.2. Lower bound evaluation. Assumption 15.21 is still used. Here our purpose is to 
evaluate a lower bound of \\uj\f — vrllr^/. The main tools are Markov and Chebyshev 
inequalities. We recall the mean value of a real function / on finite set X for a probability 
measure /i on X as follows: 

EM) = 5^/(^)/^(a^)- 

x€X 

We also recall the variance of / for ju: 

VM) = E.if) - EM)"- 

Let A be a subset of X. Then, we have the following inequality from the dehnition of the 
total variation distance: 


\Wn - ttIItv > WnM - 7i{A)\, 

where fi{A) = ^ probability measure /i. We set q 

1 "" 

q{xi....,Xn : n) = 

n 

i=i 

We identify q with a function Q on {K/LY by Q{x) = q{x). 


(5.1) 


11 




i.e., 


(5.2) 


Proposition 5.5. EYQ) = 0. 

Proof. The orthogonality relation of the zonal spherical functions and the formula fl5.2p 
gives 

EniQ) = — ^ Q{x) = 0 . 

xe{KIL)^ 

□ 
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Proposition 5.6. Set uj = Then we have 


1 ^ 1 '\ 

Q'^ix) = — + - ^ 1 -) ^{n-2,2fln-2){.^) ■ 

irt 7% \ 71/ j 

J=l 

Here 0^ = 0,..., 0. 

j 

Proof. Some properties of the zonal spherical functions give us Oj >0 (0 < j < s — 1), 
oo + • • • + Os-i = 1 (and oq = ^). The coefficient of in ^(K,L-,n){x) is given 

by '^^=i^j{xi) = nr2(n-i,i,o,...,o)(a;)- Similarly the coefficient of tQ~H\ in is 

,^o)(a;). Together, these give the claim of the propo¬ 
sition. □ 


gp 

n 


Proposition 5. T. VM = 

Proof. Proposition 15.61 enable us to compute 

V^{q) = K(Q') - E^iQf = E^iQ^) 


xa{KILY 


2 / \ 

‘x) = —. 
n 


□ 


Theorem 5.8. Let 0 < c < log?T,(r — 1) and N = 2 mp(i-M) (iQg^(^ “ 1) “ c) . If n is 
sufficiently large, then there is a constant 5 > 0, which does not depend on n such that 

,, 

\Wn — 71 TV >1-7- 


Proof. From Theorem 14.11 the iV-step probability is given by 

s—1 




k£X{s^n) 


ki 


(1 — mp) + mp —ojAxq) 
n 


j=0 


N 




where = dQ° ■ ■ ■ d^"_f. The orthogonality relations of flfcS and Proposition 15.61 give us 

mp 


Ep.*N{Q) = 


n 


-(1 - a;i(xo)) 


N 


> 0 


and 




5 — 1 


j=i 


n L 


mp 


N 


1 - —{1 - UJj{xo)) +( 1 --) 
n J n 


n 


1 N 


The following inequalities hold for n > 4: 
V,^{Q) = E,^{Q^) - 
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s—1 


1 N 


E- l-^(l--.(xo)) +(1-- 

n I n J n 

j=0 

N 1 

n 




n N 


n 


1 _ _ wi(xo)) 

n J 


^ ^ 1 ~ Qq 

“ n n 


n 


1 _ _ a;i(a;o)) 

n 


1 2Af 


On 1 — On 1 

< — H-^ < -• 

n n n 


Here we set, for 0 < c < logn(r — 1), 


N = 


n 


2mp{l — M) 


(log?T.(r — 1) — c). 


Then, we have 


mp 


E,^{Q)= ^l-^(l-a;i(xo)) 


N 


= exp 

= exp 


log (l - ^(1 - M)) N 


n 


P ^ ^ _i ^ 

= [n[r — 1)J 2 e 2 exp 
> [n[r — 1)J 262 exp 


/mp 

A 4n 


(1 — M){c — logn(r — 1))) 6 — M) 




n 


— (1 — M)(c — logn(r — 1)) ) 6 ( —(1 — M) 


8n 


2n 


Here, e{x) = jg ^ positive and monotonically increasing fnnction on (0,1) 

with hm 2 ._j.oe(x) = 1. Therefore, the exp-section of the last formnla is a monotonically 
increasing fnnction for snfhciently large n. Therefore, there is 0 < 7 < 1 for snfficiently 
large n snch that 


1 

EuAQ) > 'y[n{r - l)]" 2 e— 

Set = {x G {K/LY I ^ }• From Markov’s ineqnality, we have 

\/(r-l)n 

<A,) > 1 - ^£.( 7 ) = 1 - 

The triangle ineqnality gives us C {x G {K/LY \ |Q(a:^) — Eyj^[Q)\ > \Ey^{Q)\ — 
. ^ }■ Setting (3 = 7 e'’/^/2, we have 

-^(r-l)n 

/ A \ ^ ~ 1) 

from Chebyshev’s inequality . Therefore, from fl5.ll) . we have 


I II 1 

ll/N — TT TV > 1 “ ^ 

13^ 


jr-l) 


= 1 


4(ao + (r - 1 )) 


72 gC 


□ 


2N 







































12 


MIZUKAWA 


We now have upper and lower bounds for | li/^r ~ 7r| |tv- From the dehnition of the cut-off 
([II Dehnition 2.5.1]), we obtain the following theorem. 

Theorem 5.9. Fix a Gelfand pair {K, L) with the zonal spherical functions that are real¬ 
valued functions. Under Assumption 15.M the stochastic models defined by Deftnition \3.1\ 
have the cut-off. 

Remark 5.10. Assumption \5.A -f3) is needed for technical reasons in the paper. However, 
even for the mp >1/2 case, similar results are expected to hold. For concrete examples 
that are considered in [2] and laiz], the authors show the cut-off by taking mp = n/{n + 
1 ) > 1 / 2 . 
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